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ABSTRACT 


Properties  of  the  phase  space  of  the  standard  map  with  memory  are  investigated  This  map  was  obtained 
from  a  kicked  fractional  differential  equation.  Depending  on  the  value  of  the  map  parameter  and  the 
fractional  order  of  the  derivative  in  the  original  differential  equation,  this  nonlinear  dynamical  system 
demonstrates  attractors  (fixed  points,  stable  periodic  trajectories,  slow  converging  and  slow  diverging 
trajectories,  ballistic  trajectories,  and  fractal-like  structures)  and/or  chaotic  trajectories.  At  least  one  type 
of  fractal-like  sticky  attractors  in  the  chaotic  sea  was  observed. 

©  2009  Elsevier  B  V.  All  rights  reserved. 


1.  Introduction 

The  standard  map  (SM)  can  be  derived  from  the  differential 
equation  describing  kicked  rotator.  The  description  of  many  physi¬ 
cal  systems  and  effects  (Fermi  acceleration,  comet  dynamics,  etc.) 
can  be  reduced  to  the  studying  of  the  SM  [1],  The  SM  provides  the 
simplest  model  of  the  universal  generic  area  preserving  map  and  it 
is  one  of  the  most  widely  studied  maps.  The  topics  examined  in¬ 
clude  fixed  points,  elementary  structures  of  islands  and  a  chaotic 
sea,  and  fractional  kinetics  (1-3). 

It  was  recently  realized  that  many  physical  systems,  includ¬ 
ing  systems  of  oscillators  with  long  range  interaction  [4,5],  non- 
Markovian  systems  with  memory  ([6,  Chapter  10),  [7-11]),  fractal 
media  |12),  etc.,  can  be  described  by  the  fractional  differential 
equations  (FDE)  (6,13,14].  As  with  the  usual  differential  equations, 
the  reduction  of  FDEs  to  the  corresponding  maps  can  provide  a 
valuable  tool  for  the  analysis  of  the  properties  of  the  original  sys¬ 
tems.  As  in  the  case  of  the  SM,  the  fractional  standard  map  (FSM), 
derived  in  |15]  from  the  fractional  differential  equation  describing 
a  kicked  system,  is  perhaps  the  best  candidate  to  start  a  general  in¬ 
vestigation  of  the  properties  of  maps  which  can  be  obtained  from 
FDEs. 


*  Corresponding  author  at:  Courant  Institute  of  Mathematical  Sciences,  New  York 
University,  2S1  Mercer  St.,  New  York.  NY  10012,  USA. 
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As  it  was  shown  in  ]  1 5],  maps  that  can  be  derived  from  FDEs 
are  of  the  type  of  discrete  maps  with  memory.  One-dimensional 
maps  with  memory,  in  which  the  present  state  of  evolution  de¬ 
pends  on  all  past  states,  were  studied  previously  in  (16-21).  They 
were  not  derived  from  differential  equations.  Most  results  were  ob¬ 
tained  for  the  generalizations  of  the  logistic  map. 

In  the  physical  systems  the  transition  from  integer  order  time 
derivatives  to  fractional  (of  a  lesser  order)  introduces  additional 
damping  and  is  similar  in  appearance  to  additional  friction  ]6,22|. 
Accordingly,  in  the  case  of  the  FSM  we  may  expect  transformation 
of  the  islands  of  stability  and  the  accelerator  mode  islands  into  at¬ 
tractors  (points,  attracting  trajectories,  strange  attractors).  Because 
the  damping  in  systems  with  fractional  derivatives  is  based  on  the 
internal  causes  different  from  the  external  forces  of  friction  (22, 
23],  the  corresponding  attractors  are  also  different  from  the  attrac¬ 
tors  of  the  regular  systems  with  friction  and  are  called  fractional 
attractors  [22].  Even  in  one-dimensional  cases  (16-21]  most  of  the 
results  were  obtained  numerically.  An  additional  dimension  makes 
the  problem  even  more  complex  and  most  of  the  results  in  the 
present  Letter  were  obtained  numerically. 

2.  FSM,  initial  conditions 

The  standard  map  in  the  form 
Pn+1  =Pn  “  KsinXn, 

Xn  +  1=Xn  +  Pn  +  1  (mod27T)  (1) 
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can  be  derived  from  the  differential  equation 

x+Ksin(x)]T«5^— n)=0.  (2) 

n=0  '  ' 

By  replacing  the  second-order  time  derivative  in  Eq.  (2)  with 
the  Riemann-Liouville  derivative  one  obtains  a  fractional 
equation  of  motion  in  the  form 

0D?x  +  Ksin(x)£<$(L-n)=0  (l<as$2),  (3) 

n=0  '  ' 


where 


oD?x(t)  =  D?01?-ax(t) 

1  dm  f  x(r)dr 

~  r(m  —  a)  dt «  J  (t  -  r)a-m+1 

o 


(m  -  1  <  a  <  m), 
(4) 


D™  ~dm/dtm,  and  o/"  is  a  fractional  integral.  The  initial  conditions 
for  (3)  are 


(oD“  ’x)(0+)=Pi. 

(0D?-2x)(0+)=b.  (5) 

The  Cauchy  type  problem  (3)  and  (5)  is  equivalent  to  the  Volterra 
integral  equation  of  the  second  kind  [24-26] 


x(0  = 


Pi  fg-l  ,  ^  ta-2 

T(a)  r(a-l) 

K  r  sin[x(r)]  ~n)dr 

r (a)J  (t  -  r)]~a 


0 

Defining  the  momentum  as 


(6) 


p(t)  =  0D“-1x(t),  (7) 

and  performing  integration  in  (6)  one  can  derive  the  equation  for 
the  FSM  in  the  form  (for  the  thorough  derivation  see  [26]) 


Pn+l  =Pn  -Ksinxn,  (8) 

1  " 

Xn+1  =  —  ]Tpi+,Va(n-i  +  l) 

V  '  1=0 

+  ■  -  '  .  (n  +  l)g~2  (mod 271-),  (9) 

T  (or  —  1 ) 

where 

Va(m)  =  m“-1  —  (m  —  1)“_1.  (10) 

Here  it  is  assumed  that  T  =  1  and  1  <  a  <  2.  The  form  of  Eq.  (9) 
which  provides  a  more  clear  correspondence  with  the  SM  (a  =2) 
in  the  case  b  =  0  is  presented  in  Section  4  (Eq.  (31)). 

The  second  initial  condition  in  (5)  can  be  written  as 


in  the  limiting  case  a  =  2  coincide  with  the  equations  for  the  stan¬ 
dard  map  under  the  condition  xo  =0. 

In  this  Letter  the  FSM  is  taken  in  the  form  derived  in  [15]  which 
coincides  with  (8)  and  (9)  if  b  =  0.  It  is  also  assumed  that  Xo  =  0 
and  the  results  can  be  compared  to  those  obtained  for  the  SM  with 
xo  =  0  and  arbitrary  po-  As  a  test,  for  the  SM  and  for  the  FSM  with 
a  =  2  and  the  same  initial  conditions  numerical  calculations  show 
that  phase  portraits  look  identical. 

System  of  Eqs.  (8)  and  (9)  can  be  considered  either  in  a  cylin¬ 
drical  phase  space  (x  mod  27t)  or  in  unbounded  phase  space.  The 
second  case  is  convenient  to  study  transport.  The  trajectories  in 
the  second  case  are  easily  related  to  the  first  case.  The  FSM  has 
no  periodicity  in  p  (the  SM  does)  and  cannot  be  considered  on  a 
torus. 

3.  Stable  fixed  point 

The  SM  has  stable  fixed  points  at  (0, 27rn)  for  K  <  Kc  =  4.  It  is 
easy  to  see  that  point  (0,0)  is  also  a  fixed  point  for  the  FSM.  Di¬ 
rect  computations  using  (8)  and  (9)  demonstrate  that  for  the  small 
initial  values  of  po  there  is  a  clear  transition  from  the  convergence 
to  the  fixed  point  to  divergence  when  the  value  of  the  parameter 
K  crosses  the  curve  K  —  Kc(a )  on  Fig.  1(a)  from  smaller  to  larger 
values. 

The  following  system  describes  the  evolution  of  trajectories 
near  fixed  point  (0, 0) 


<$Pn+l  —  &Pn  ~  K&Xn, 


1  n ' 

<$*n+1  =  —  i  +  !)• 


i=0 

The  solution  can  be  found  in  the  form 

<n>« 


<5x„  = 


Po 

Ha) 


VaiK_\ 

ana)) 


(n  >  0). 


(12) 

(13) 

(14) 

(15) 


The  origin  of  the  terms  in  parentheses,  as  well  as  the  definition 

oo 


k=  1 


(16) 


will  become  clear  in  Section  5.  Eqs.  (12)-(16)  lead  to  the  following 
iterative  relationships 

n 

xn+i,i  =  -53(n-m  +  l)"_1xm.i_,  (0  <  i  <  n),  (17) 

m=i 
n 

P  n+i,i  =  -]Txmf_i  (0<t<n)  (18) 


with  the  initial  and  boundary  conditions 


(0D“_2x)(0+)  = 


lim  o J2  “x(f) 
r->0+ 


=  lim  - 

t— ►0+  P  (2 


-ct)  J 


x(r)dr 
(t  -  r)°'— ’ 


=  b  (1  <  or  <  2), 


(11) 


which  requires  b  =  0  in  order  to  have  a  solution  bounded  at  t  =  0 
for  a  <2.  The  assumption  b  =  0  leads  to  the  FSM  equations  which 


Pn+I.n  —  Xn+1,n  —  (— l)n>  Pn+1,0  —  1> 

Xn+I.0  =  (n+1)°'-1.  (19) 

From  (17)  and  (18)  it  is  clear  that  the  series  (14)  and  (15)  are 
alternating  and  it  is  natural  to  apply  the  Dirichlet’s  test  to  verify 
their  convergence.  This  can  be  done  by  considering  the  totals 


(20) 
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Fig.  1.  Stability  of  the  fixed  point  (0,0):  (a)  The  fixed  point  is  stable  below  the 
curve  K  =  Kc(a)\  (b)  Values  of  S<x  and  l0 0  obtained  after  20000  iterations  of 
Eqs.  (22)  and  (23).  As  a  -*•  2  the  values  S and  /0 0  increase  rapidly.  For  a  =  1 .999, 
Soo  w  276  and  %  552  after  20000  iterations;  (c)  An  example  of  the  typical  evo¬ 
lution  of  Soo  and  /<x>  over  the  first  200  iterations  for  1  <  a  <  2  This  particular 
figure  corresponds  to  a  =  1.8;  (d)  Deviation  of  the  values  Sn  and  /„  from  the  val¬ 
ues  Soo  **  2.04337  and  l0 0  «  3.37416  for  a  =  1.8  during  the  first  20000  iterations 
(this  type  of  behavior  remains  for  1  <  a  <  2);  (e)  Evolution  of  trajectories  with 
Po  =  1.5  +  0  0005i,  0  ^  i  <  200  for  the  case  K  =  3,  a  =  1 .9.  The  line  segments  cor¬ 
respond  to  the  nth  iteration  on  the  set  of  trajectories  with  close  initial  conditions. 
The  evolution  of  the  trajectories  with  smaller  p0  is  similar;  (0  10s  iterations  on 
both  of  two  trajectories  for  K  -  2,  a  =  1.4.  The  one  at  the  bottom  with  po  =  0,3 
is  a  fast  converging  trajectory  The  upper  trajectory  with  p0  =  5.3  is  an  example  of 
the  ASCT  in  which  Piooooo  **  0.042. 


They  obey  the  following  iterative  rules 


n— 1 


Sn  =  n“  1  -  —  Yj(n  -  i)“  1  Si,  S,  =  l, 

*  al  . 


i=1 


(22) 

(23) 


Computer  simulations  show  that  values  of  Sn  and  /n  converge 
to  the  values  (-l)n+1Soo  and  (-l)n+l/oo  depicted  on  Fig.  1(b). 
Figs.  1(c),  (d)  show  an  example  of  the  typical  evolution  of  Sn  and 
ln  over  the  first  20000  iterations.  It  means  that  the  condition  of 
convergence  of  Spn  and  Sxn  is 


£si<i 


(24) 


Numerical  evaluation  of  the  equality  K  =2T (a)/Vai  perfectly  re¬ 
produces  the  curve  on  Fig.  1(a)  obtained  by  the  direct  computa¬ 
tions  of  (8)  and  (9). 

Because  not  only  the  stability  problem  (12)  and  (13),  but  also 
the  original  map  (8)  and  (9),  contains  convolutions,  the  use  of  gen¬ 
erating  functions  [27],  which  allows  transformations  of  sums  of 
products  into  products  of  sums,  could  be  utilized  in  the  inves¬ 
tigation  of  the  FSM  and  some  other  maps  with  memory.  As  an 
example,  in  the  particular  case  of  the  stability  problem  (12)  and 
(13),  the  introduction  of  the  generating  functions 


K  00 

w«(o=— -y>+ D"-1- 

r(a)  ^ 

00 

X(t)  = 

i=0 

00 

P(t)  =  ^Spiti, 

1=0 

leads  to 


(25) 

(26) 
(27) 


X(t)  = 


PoWa(() 


t 


P(t)  =  Po 


K  l-t(l-Wff(t)) 
1  +  W„(t) 


1  -t(l  -wa(t)) 


(28) 

(29) 


Now  the  original  problem  is  reduced  to  the  problem  of  the  asymp¬ 
totic  behavior  at  t  =  0  of  the  derivatives  of  the  analytic  functions 
X(t)  and  P(t),  which  is  still  quite  complex  and  is  not  considered 
in  this  Letter. 

In  the  region  of  the  parameter  space  where  the  fixed  point  is 
stable,  the  fixed  point  is  surrounded  by  a  finite  basin  of  attrac¬ 
tion,  whose  width  W  depends  on  the  values  of  K  and  a.  For 
example,  for  K  =3  and  a  =  1.9  the  width  of  the  basin  of  at¬ 
traction  is  1.6  <  W  <1.7.  Simulations  of  thousands  of  trajectories 
with  po  <  1.6  performed  by  the  authors,  of  which  only  200  (with 
1.5  <  po  <  1.6)  are  presented  in  Fig.  1(e),  show  only  converging 
trajectories,  whereas  among  200  trajectories  with  1.6  <  po  <  1.7 
in  Fig.  2(a)  there  are  trajectories  converging  to  the  fixed  point  as 
well  as  some  trajectories  converging  to  attracting  slow  diverging 
trajectories  (ASDT),  whose  properties  will  be  discussed  in  the  fol¬ 
lowing  section.  Trajectories  in  Fig.  1(e)  converge  very  rapidly,  ln 
the  case  K  =  2  and  a  =  1.4  in  addition  to  the  trajectories  which 
converge  rapidly  and  ASDTs  there  exist  attracting  slow  converging 
trajectories  (ASCT)  (Fig.  1(f)). 


4.  Attracting  slow  diverging  trajectories  (ASDT) 


As  it  can  be  seen  from  Fig.  2(a),  the  phase  portrait  on  a  cylinder 
of  the  FSM  with  K  =  3  and  a  =  1  9  contains  only  one  fixed  point 
and  ASDTs  approximately  equally  spaced  along  the  p-axis.  This  re¬ 
sult  corresponds  to  the  fact  that  the  standard  map  with  K  =  3 
has  only  one  central  island.  More  complex  structure  of  the  stan¬ 
dard  map’s  phase  space  for  smaller  values  of  K  (for  example  for 
K  =  2  and  K  =  0.6)  can  explain  more  complex  structure  of  the 
FSM’s  phase  space,  where  periodic  attracting  trajectories  with  pe¬ 
riods  T  =  4  (Fig.  2(b)),  T  =  2,  and  7  =  3  (Fig.  2(c))  are  present. 

Each  ASDT  has  its  own  basin  of  attraction  (see  Fig.  2(d)).  Be¬ 
tween  those  basins  two  initially  close  trajectories  at  first  diverge, 
but  then  converge  to  the  same  or  different  fixed  point  or  ASDT. 

Numerical  evaluation  shows  that  for  ASDTs  which  converge  to 
trajectories  along  the  p-axis  (x  -*•  Xnm  =  0)  in  the  area  of  stabil¬ 
ity  (which  is  the  same  as  for  the  stability  of  the  fixed  point)  the 
following  holds  (for  large  n  see  Fig.  3(a)) 


2T  (a) 
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Fig.  2.  Phase  space  with  ASDTs:  (a)  The  same  values  of  parameters  as  in  Fig.  1(e)  but  po  =  1  6  +  0  0005i,  (b)  200  iterations  on  trajectories  with  po  =  4  +  0.02i,  0  ^  i  <  500 
for  the  case  K  =2,  or  =  1.9.  Trajectories  converging  to  the  fixed  point,  ASDTs  with  x  =  0,  and  period  4  attracting  trajectories  are  present;  (c)  2000  iterations  on  trajectories 
with  po  —  2  +  0.04i,  0  ^  i  <  50  for  the  case  K  =  0.6,  a  =  1.9.  Trajectories  converging  to  the  fixed  point,  period  2  and  3  attracting  trajectories  are  present;  (d)  The  same 
values  of  parameters  as  in  Fig.  1(e)  but  po  =  5  4-  0.005i. 


Fig.  3.  Evaluation  of  the  behavior  of  the  ASDTs:  (a)  Momenta  for  two  ASDTs  with 
xn  «==  2tt n  in  the  unbounded  space  (in  this  example  K  =  2).  The  solid  line  is  related 
to  a  trajectory  with  or  =  1  9  and  its  slope  is  01.  The  dashed  line  corresponds  to 
a  trajectory  with  a  =  1.5  and  its  slope  is  0.5;  (b)  Deviation  of  momenta  from  the 
asymptotic  formula  for  two  ASDTs  with  xn  %  2tc n  in  the  unbounded  space,  or  =  1.9, 
and  K  —  2.  The  dashed  line  has  po  =  7  and  the  solid  one  po  =  6;  (c)  Relative  de¬ 
viation  of  the  momenta  for  the  trajectories  in  (b)  from  the  asymptotic  formula; 
(d)  Deviation  of  the  x-coordinates  for  the  trajectories  in  (b)  from  the  asymptotic 
formula. 


constant  step  in  x  in  the  unbounded  space.  Then  Eq.  (9)  with  b  =  0 
gives 


1  n 

xn+l  -  X„  =  ——  J2(Pk+ 1  -  Pk) Va(n  -  k  +  1) 

r(a)  t[ 


+ 


Pi 

r(a) 


v*(n  +  l). 


(31) 


For  large  n  the  last  term  is  small  na~2)  and  the  following  holds 


n 

1  -  Pk) Va(n  -  k  +  1)  =  2nMr(a).  (32) 

1 

With  the  assumption  pn  —  n2~a  it  can  be  shown  that  for  values 
of  a  >  1.8  considered  the  terms  in  the  last  sum  with  large  k  are 
small  and  in  the  series  representation  of  Va{n  —  k  + 1)  it  is  possible 
to  keep  only  terms  of  the  highest  order  in  k/n.  Thus,  (32)  leads  to 
the  approximations 


2nMr{ct)n2~a 

Pn^P  o  H - ; - 

a  - 1 

27rM(2  -of)r(of) 

Xn*"  X(a  —  l)n“-' 


(33) 

(34) 


In  the  case  K  =  2,  of  =  1.9  Figs.  3(b)-(d)  show  two  trajectories  with 
M  =  1  (initial  momenta  po  =  6  and  po  =  7)  approaching  an  ASDT: 
the  deviation  from  the  asymptotic  (33)  and  (34)  and  the  relative 
difference  with  respect  to  (33). 


5.  Period  2  stable  trajectory 


pn  =  Cn2-a.  (30)  The  SM  has  two  stable  points  of  the  period  7  =  2  trajectory  for 

4  <  K  <  2tt  with  the  property 

The  constant  C  can  be  easily  evaluated  for  1.8  <  a  <  2.  Consider 
an  ASDT  with  X/im  =  0,  7  =  1,  and  2 ttM,  where  M  is  an  integer,  pn+ 1  =  —  pn, 


xn+i  =  -xn. 


(35) 
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Fig.  4.  Period  2  stable  trajectory:  (a)  An  example  of  T  =  2  attractor  for  K  =  4.5, 
a  =  1.9.  One  trajectory  with  xq  =  0,  po  =  0.513;  (b)  pi  of  x/  for  the  case  of  K  =  4  5, 
(c)  pi  of  a  for  the  case  of  K  =  4.5;  (d)  X/  of  a  for  the  case  of  K  =  4  5;  (e)  p„  -  pi 
for  the  trajectory  in  (a).  After  1000  iterations  |pn  -  pH  <  10~7;  (f)  x«  -x*  for  the 
trajectory  in  (a).  After  1000  iterations  |xn  -  xi|  <  10  7 . 

The  same  points  persist  in  the  numerical  experiments  for  the  FSM 
(Fig.  4(a)).  These  points  are  attracting  most  of  the  trajectories  with 
small  po.  Assuming  the  existence  of  a  T  =  2  attracting  trajectory, 
it  is  possible  to  calculate  the  coordinates  of  its  attracting  points 
(Xj,  pj)  and  (—x/,  — p/).  In  this  case  from  (8)  and  (9) 


Pi  =  jsin(x,),  (36) 

K 

x,  =  ——sln(xi)T(-1)k+'Va(k).  (37) 

2r(co  f~l 


Fig.  5.  Cascade  of  bifurcation  type  trajectories:  (a)  120000  iterations  on  a  single 
trajectory  with  K  =  4.5,  a  m  1.65,  po  =  0.3.  The  trajectory  occasionally  sticks  to 
a  CBTT  but  then  always  recovers  into  the  chaotic  sea;  (b)  100000  iterations  on  a 
trajectory  with  K  -  3.5,  a  =  1.1,  po  =  20.  The  trajectory  very  fast  turns  into  a  CBTT 
which  slowly  converges  to  a  fractal  type  area. 


Finally,  the  equation  for  Xj  takes  the  form 

Xl=  2f(a)ValSm{Xl)'  (38) 

where 


K  >  4  for  the  SM.  The  results  of  calculations  of  the  X/  and  pj  for 
the  cases  K  =  4.5,  1  <  a  <  2  presented  in  Figs.  4(b)— (d)  perfectly 
coincide  with  the  results  of  the  direct  computations  of  (8)  and  (9) 
with  b  =  0.  After  1000  iterations  presented  in  Figs.  4(e),  (f)  the 
values  of  deviations  |pn  —  P/1  and  |xn  —  X/ 1  are  less  than  10“7. 


oo 

Va,=  £(-l)k+1Va(J0  (39) 

k=l 

and  can  be  easily  calculated  numerically.  From  (38)  the  condition 
of  the  existence  of  T  =  2  trajectory 

ir>ir«W-^.  (40) 

Val 

is  exactly  opposite  to  (24).  It  is  satisfied  above  the  curve  K  =  Kc{a) 
on  Fig  1(a).  For  a  =  2  (40)  produces  the  well-known  condition 


6.  Cascade  of  bifurcations  type  trajectories  (CBTT) 

Period  2  stable  trajectories  have  limited  basins  of  attraction. 
Trajectories  that  don’t  fall  into  those  areas  reveal  a  diverse  vari¬ 
ety  of  properties,  from  period  two  slow  attracting  trajectories  to 
fractal  type  attractors  and  cascade  of  bifurcations  type  trajectories 
(CBTT).  Fig.  5(a)  presents  a  single  chaotic  trajectory  which  sticks 
to  the  areas  similar  to  the  cascade  of  bifurcations  which  are  well 
known  for  the  logistic  map.  In  Fig.  5(b)  a  single  trajectory  falls  very 
rapidly  into  one  of  the  attracting  CBTTs.  Because  the  bifurcation  di- 
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Fig.  6.  Examples  of  phase  space  for  K  >  2n\  (a)  An  attracting  ballistic  trajectory  with 

agram  of  the  logistic  map  has  fractal  properties  (see  for  example 
Chapter  2  in  [28]),  it  is  expected  that  the  structure  to  which  this 
trajectory  slowly  converges  also  possesses  fractal  features. 

The  properties  of  this  type  of  attractors,  as  well  as  the  prop¬ 
erties  different  types  of  observed  during  computer  simulations 
chaotic,  attracting,  and  ballistic  trajectories  for  K  >  2n  (see  Fig.  6) 
will  be  considered  in  the  subsequent  article. 

7.  Fractional  attractors  and  their  stability 

The  problems  of  existence  and  stability  of  the  fractional  attrac¬ 
tors  for  the  systems  described  by  the  FDEs  were  addressed  in  a 
few  recent  papers.  It  was  noticed  in  [22]  that  the  properties  of  the 
fractional  chaotic  attractors  are  different  from  the  properties  of  the 
"regular"  chaotic  attractors  and  may  have  some  pseudochaotic  fea¬ 
tures.  The  problem  of  existence  of  multi-scroll  fractional  chaotic 
attractors  was  considered  in  [29].  The  problem  of  stability  of  the 
stationary  solutions  (fixed  points  for  ODEs)  of  systems  described 
by  the  fractional  ODEs  and  PDEs  was  considered  in  [30-32].  In 
the  above  mentioned  articles  the  equations  contained  the  Caputo 
fractional  derivatives,  whereas  in  the  present  Letter  the  Riemann- 
Liouville  fractional  derivative  is  used.  This  fact  does  not  allow  a 
direct  comparison  of  the  results.  The  results  [22,29-32]  were  sup¬ 
ported  by  a  relatively  small  number  of  computations  and  this  is 
understandable,  taking  into  account  all  the  difficulties  of  perform¬ 
ing  numerical  simulations  for  the  equations  with  fractional  deriva¬ 
tives. 

The  use  of  the  FSM,  which  is  equivalent  to  the  original  FDE, 
allows  performing  thousands  of  runs  of  simulations  of  the  kicked 
fractional  system  with  two  parameters:  K  and  a.  The  FSM  also 
allows  making  some  analytic  deductions  and  revealing  some  prop¬ 
erties  of  the  fractional  attractors  which  were  not  reported  before: 

(a)  The  stability  of  the  fixed  point  (0,0)  of  the  FSM  is  different 
not  only  from  the  stability  of  the  fixed  point  in  the  domain  of  the 
regular  motion  (zero  Lyapunov  exponent)  of  the  SM,  but  also  from 
the  stability  of  fixed  attracting  points  of  the  regular  (not  fractional) 
dissipative  systems  like,  for  example,  the  dissipative  standard  map 
(Zaslavsky  map)  ]33].  The  difference  is  in  the  way  in  which  trajec¬ 
tories  approach  the  attracting  point.  In  the  FSM  this  way  depends 
on  the  initial  conditions.  For  example,  in  Fig.  1(f)  there  are  two 
trajectories  approaching  the  same  fixed  point:  one  is  fast  spiraling 
into  the  attractor  and  the  other  is  slowly  converging. 

(b)  Stable  period  2  attracting  trajectories  exist  only  in  the 
asymptotic  sense— they  do  not  represent  any  real  periodic  solu¬ 
tions.  If  the  initial  condition  is  chosen  in  a  period  two  stable 
attracting  point,  this  trajectory  will  immediately  jump  out  of  this 
point  and  where  it  will  end  depends  on  the  values  of  K  and  a . 

(c)  All  the  FSM  attractors  exist  in  the  sense  that  there  are 
trajectories  which  converge  into  those  attractors.  But  if  an  initial 
condition  is  taken  on  any  of  the  attracting  trajectories  (except  for 


1500  - 


=  6  908745,  a  =  1.999,  po  =0.7;  (b)  A  chaotic  trajectory  for  K  =  6  908745,  a  =  1.9. 

the  fixed  point),  they  will  most  likely  not  evolve  along  the  same 
trajectory. 

8.  Conclusion 

In  this  Letter  properties  of  the  phase  space  of  the  FSM  were 
investigated.  It  was  shown  that  islands  of  regular  motion  of  the 
SM  in  the  FSM  turn  into  attractors  (points,  attracting  trajecto¬ 
ries,  and  fractal-like  structures).  Properties  of  the  attracting  fixed 
points,  period  two  trajectories,  ASCTs,  and  ASDTs  were  considered. 
This  consideration  allows  the  description  of  the  evolution  of  the 
dynamical  variable  x  of  the  original  fractional  dynamical  system, 
a  system  described  by  the  FDE  reducible  to  the  FSM. 

The  explanation  of  the  CBTTs,  which  are  interesting  phenom¬ 
ena,  requires  further  detailed  investigation.  Chaotic  trajectories  that 
spend  some  time  near  CBTTs,  which  can  be  called  "sticky  attrac¬ 
tors"  in  analogy  to  "sticky  islands"  of  the  SM,  are  good  candidates 
for  the  investigation  of  anomalous  diffusion.  Phase  space  transport 
was  not  considered  in  this  Letter.  How  general  the  properties  of 
the  phase  space  of  the  FSM  are  will  become  clear  after  further  in¬ 
vestigations  of  different  fractional  maps,  maps  with  memory  which 
can  be  derived  from  the  FDEs,  and  particular  those  suggested  in 
[15],  will  be  conducted.  The  fact  that  so  many  physical  systems  can 
be  reduced  to  studying  of  the  SM  gives  a  hope  that  those  physi¬ 
cal  systems  which  can  be  reduced  to  studying  the  FSM  will  be 
found. 
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Abstract 

Discrete  maps  with  long-term  memory  are  obtained  from  nonlinear  differential 
equations  with  Riemann-Liouville  and  Caputo  fractional  derivatives.  These 
maps  are  generalizations  of  the  well-known  universal  map.  The  memory  means 
that  their  present  state  is  determined  by  all  past  states  with  special  forms  of 
weights.  To  obtain  discrete  maps  from  fractional  differential  equations,  we 
use  the  equivalence  of  the  Cauchy-type  problems  and  to  the  nonlinear  Volterra 
integral  equations  of  the  second  kind.  General  forms  of  the  universal  maps 
with  memory,  which  take  into  account  general  initial  conditions  for  the  cases 
of  the  Riemann-Liouville  and  Caputo  fractional  derivative,  are  suggested. 

PACS  numbers:  05.45. -a,  45.10. Hj 
Mathematics  Subject  Classification:  26A33,  37E05 


1.  Introduction 

A  dynamical  system  consists  of  a  set  of  possible  states,  together  with  a  rule  that  determines  the 
present  state  in  terms  of  past  states.  If  we  require  that  the  rule  be  deterministic,  then  we  can 
define  the  present  state  uniquely  from  the  past  states.  A  discrete-time  system  without  memory 
takes  the  current  state  as  input  and  updates  the  situation  by  producing  a  new  state  as  output. 
All  physical  classical  models  are  described  by  differential  or  integro-differential  equations, 
and  the  discrete-time  systems  can  be  considered  as  a  simplified  version  of  these  equations.  A 
discrete  form  of  the  time  evolution  equation  is  called  the  map.  Maps  are  important  because 
they  encode  the  behavior  of  deterministic  systems.  The  assumption  of  determinism  is  that  the 
output  of  the  map  can  be  uniquely  determined  from  the  input.  In  general,  the  present  state  is 
uniquely  determined  by  all  past  states,  and  we  have  a  discrete  map  with  memory.  Discrete 
maps  are  used  for  the  study  of  evolution  problems,  possibly  as  a  substitute  of  differential 
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equations  [1-3].  They  lead  to  a  much  simpler  formalism,  which  is  particularly  useful  in 
simulations.  The  universal  discrete  map  is  one  of  the  most  widely  studied  maps.  In  this  paper, 
we  consider  discrete  maps  with  memory  that  can  be  used  to  study  solutions  of  fractional 
differential  equations  [4-7], 

The  nonlinear  dynamics  can  be  considered  in  terms  of  discrete  maps.  It  is  a  very  important 
step  in  understanding  the  qualitative  behavior  of  systems  described  by  differential  equations. 
The  derivatives  of  non-integer  orders  are  a  generalization  of  the  ordinary  differentiation  of 
integer  order.  Fractional  differentiation  with  respect  to  time  is  characterized  by  long-term 
memory  effects.  The  discrete  maps  with  memory  are  considered  in  [8-14].  The  interesting 
question  is  a  connection  of  fractional  differential  equations  and  discrete  maps  with  memory. 
It  is  important  to  derive  discrete  maps  with  memory  from  the  equation  of  motion. 

In  [14],  we  prove  that  the  discrete  maps  with  memory  can  be  obtained  from  differential 
equations  with  fractional  derivatives.  The  fractional  generalization  of  the  universal  map  was 
obtained  [14]  from  a  differential  equation  with  Riemann-Liouville  fractional  derivatives.  The 
Riemann-Liouville  derivative  has  some  notable  disadvantages  such  as  the  hyper-singular 
improper  integral,  where  the  order  of  singularity  is  higher  than  the  dimension,  and  nonzero  of 
the  fractional  derivative  of  constants,  which  would  entail  that  dissipation  does  not  vanish  for  a 
system  in  equilibrium.  The  desire  to  formulate  initial  value  problems  for  mechanical  systems 
leads  to  the  use  of  Caputo  fractional  derivatives  rather  than  the  Riemann-Liouville  fractional 
derivative. 

It  is  possible  to  state  that  the  Caputo  fractional  derivatives  allow  us  to  give  more 
clear  mechanical  interpretation.  At  the  same  time,  we  cannot  state  that  the  Riemann- 
Liouville  fractional  derivative  does  not  have  a  physical  interpretation  and  that  it  shows 
unphysical  behavior.  Physical  interpretations  of  the  Riemann-Liouville  fractional  derivatives 
are  more  complicated  than  Caputo  fractional  derivatives  But  the  Riemann-Liouville  fractional 
derivatives  naturally  appear  for  real  physical  systems  in  electrodynamics.  We  note  that  the 
dielectric  susceptibility  of  a  wide  class  of  dielectric  materials  follows,  over  extended  frequency 
ranges,  a  fractional  power-law  frequency  dependence  that  is  called  the  ‘universal’  response 
[15,  16].  As  was  proved  in  [17,  18],  the  electromagnetic  fields  in  such  dielectric  media  are 
described  by  differential  equations  with  Riemann-Liouville  fractional  time  derivatives.  These 
fractional  equations  for  ‘universal’  electromagnetic  waves  in  dielectric  media  are  common  to 
a  wide  class  of  materials,  regardless  of  the  type  of  physical  structure,  chemical  composition, 
or  of  the  nature  of  the  polarizing  species.  Therefore,  we  cannot  state  that  Riemann-Liouville 
fractional  time  derivatives  do  not  have  a  physical  interpretation.  The  physical  interpretation 
of  these  derivatives  in  electrodynamics  is  connected  with  the  frequency  dependence  of  the 
dielectric  susceptibility.  As  a  result,  the  discrete  maps  with  memory  that  are  connected  with 
differential  equations  with  Riemann-Liouville  fractional  derivatives  are  very  important  to 
physical  applications,  and  these  derivatives  naturally  appear  for  real  physical  systems. 

For  computer  simulation  and  physical  application,  it  is  very  important  to  take  into  account 
the  initial  conditions  for  discrete  maps  with  memory  that  are  obtained  from  differential 
equations  with  Riemann-Liouville  fractional  time  derivatives.  In  [14],  these  conditions  are 
not  considered.  In  this  paper  to  take  into  account  the  initial  condition,  we  use  the  equivalence 
of  the  differential  equation  with  Riemann-Liouville  and  Caputo  fractional  derivatives  and  the 
Volterra  integral  equation.  This  approach  is  more  general  than  the  auxiliary  variable  method 
that  is  used  in  [14].  The  proof  of  the  result  for  Riemann-Liouville  fractional  derivatives 
is  more  complicated  in  comparison  with  the  results  for  the  Caputo  fractional  derivative.  In 
this  paper,  we  prove  that  the  discrete  maps  with  memory  can  be  obtained  from  differential 
equations  with  the  Caputo  fractional  derivative.  The  fractional  generalization  of  the  universal 
map  is  obtained  from  a  fractional  differential  equation  with  Caputo  derivatives. 
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The  universal  maps  with  memory  are  obtained  by  using  the  equivalence  of  the  fractional 
differential  equation  and  the  Volterra  integral  equation.  We  reduce  the  Cauchy-type  problem 
for  the  differential  equations  with  the  Caputo  and  Riemann-Liouville  fractional  derivatives 
to  nonlinear  Volterra  integral  equations  of  second  kind.  The  equivalence  of  this  Cauchy- 
type  problem  for  the  fractional  equations  with  the  Caputo  derivative  and  the  correspondent 
Volterra  integral  equation  was  proved  by  Kilbas  and  Marzan  in  [19,  20].  We  also  use  that 
the  Cauchy-type  problem  for  the  differential  equations  with  the  Riemann-Liouville  fractional 
derivative  can  be  reduced  to  a  Volterra  integral  equation.  The  equivalence  of  this  Cauchy-type 
problem  and  the  correspondent  Volterra  equation  was  proved  by  Kilbas,  Bonilla  and  Trujillo 
in  [21,22]. 

In  section  2,  differential  equations  with  integer  derivative  and  universal  maps  without 
memory  are  considered  to  fix  notations  and  provide  convenient  references.  In  section  3, 
fractional  differential  equations  with  the  Riemann-Liouville  derivative  and  universal  maps 
with  memory  are  discussed.  In  section  4,  the  difference  between  the  Caputo  and  Riemann- 
Liouville  fractional  derivatives  is  discussed.  In  section  5,  fractional  differential  equations 
with  the  Caputo  derivative  and  correspondent  discrete  maps  with  memory  are  considered.  A 
fractional  generalization  of  the  universal  map  is  obtained  from  kicked  differential  equations 
with  the  Caputo  fractional  derivative  of  order  1  <  a  ^  2.  The  usual  universal  map  is  a  special 
case  of  the  universal  map  with  memory.  Finally,  a  short  conclusion  is  given  in  section  6. 

2.  Integer  derivative  and  universal  map  without  memory 

In  this  section,  differential  equations  with  derivative  of  integer  order  and  the  universal  map 
without  memory  are  considered  to  fix  notations  and  provide  convenient  references. 

Let  us  consider  the  equation  of  motion 


(1) 


in  which  perturbation  is  a  periodic  sequence  of  delta-function-type  pulses  (kicks)  following 
with  period  T  =  lixjv,  K  is  an  amplitude  of  the  pulses,  D2  =  d2/dt2,  and  G[x]  is  some 
real-valued  function.  It  is  well  known  that  this  differential  equation  can  be  represented  in  the 
form  of  the  discrete  map 


(2) 


**+i  -xn  =  pn+\T,  pn+\  —  pn  =  —KT  G[xn], 


Equations  (2)  are  called  the  universal  map.  For  details,  see  for  example  [1-3]. 

The  traditional  method  of  derivation  of  the  universal  map  equations  from  the  differential 
equations  is  considered  in  section  5.1  of  [2]  We  use  another  method  of  derivation  of  these 
equations  to  fix  notations  and  provide  convenient  references.  We  obtain  the  universal  map  by 
using  the  equivalence  of  the  differential  equation  and  the  Volterra  integral  equation. 

Proposition  1.  The  Cauchy-type  problem  for  the  differential  equations 


D)x{t)  =  p(t), 


(3) 


(4) 


with  the  initial  conditions 


x(0)=;to,  p(0)  =  po 


(5) 
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is  equivalent  to  the  universal  map  equations  of  the  form 

n 

*„+i  =  *o  +  Po(n  +  1  )T-KT2J2  G[xk]  (n  +  1  -  *),  (6) 

*=1 

n 

Pn+I  =  Po  ~  KT  ^G[x*].  (7) 

k= 1 


(8) 


k=\ 

The  solution  of  the  left  side  of  the  {n  +  l)th  kick 

xn+ 1  =  x(tn+ 1  -  0)  =  lim  x(T(n  +  1)  -  e ), 

f  — ►(>+ 

Pn+1  =  P0n+\  -  0)  =  lim  p(7(n  +  1)  -£), 


(10) 


Proof.  Consider  the  nonlinear  differential  equation  of  second  order 
D}x(t)  =  G[i,x(t)],  (0  </<//) 

on  a  finite  interval  [0,  t /]  of  the  real  axis,  with  the  initial  conditions 

x(0)  =  x0,  {D'x)(  0)  =  p0,  (9) 

The  Cauchy-type  problem  of  the  form  (8),  (9)  is  equivalent  to  the  Volterra  integral  equation 
of  second  kind 

x(t)  =  Xo  +  pot  +  f  dr  G[r,  x(z)]  ( t  -  r). 

J  o 

Using  the  function 

•*(/)]  =  -KG[x(/)i: 
for  nT  <  t  <  (n  +  \)T,  we  obtain 

n 

x(t)  =  xo  +  pot  -  KT  Y, G[x(kT) ]  ( t  -  kT). 

k=\ 

For  the  momentum  p(t)  =  D/jc(0*  equation  (11)  gives 

n 

p(t)  =  Po-  ktJ2g^( kT)l 


(11) 

(12) 

(13) 

(14) 


where  tn+ \  =  (n  +  1)7,  has  the  form  (6)  and  (7). 

This  ends  the  proof.  □ 

Remark  1.  We  note  that  equations  (6)  and  (7)  can  be  rewritten  in  the  form  (2).  Using 
equations  (6)  and  (7),  the  differences  xn+\  —  xn  and  pn+\  —  pn  give  equations  (2)  of  the 
universal  map. 

Remark  2.  If  G[x]  =  —  jc,  then  equations  (2)  give  the  Anosov-type  system 

Xn+ 1  Xn  =  Pn+\  7 ,  7n+l  Pn  =  KT Xn.  (15) 

For  G[;t]  =  sin(x),  equations  (2)  are 

Xn+1  -xn  =  pn+\T,  pn+\  -  Pn  =  -KT  sin(xn).  (16) 

This  map  is  known  as  the  standard  or  Chirikov-Taylor  map  [  1  ]. 
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3.  Riemann-Liouville  fractional  derivative  and  universal  map  with  memory 


In  this  section,  we  discuss  nonlinear  differential  equations  with  the  left-sided  Riemann- 
Liouville  fractional  derivative  o Df  defined  fora  >  0  by 


Id n  jc(r)dr 

oCMO  .  Of  . /;-*»>  =  l  ")• 

where  D*  =  dn/dtn,  ando If  is  a  fractional  integration  [4,  6,  7]. 

We  consider  the  fractional  differential  equation 


(17) 


oD?x(t)  =  G[t,x(t)}, 


(18) 


where  G[r,  jc (r)]  is  a  real-valued  function,  0  ^  n  —  1  <  a  ^  n,  and  t  >  0,  with  the  initial 
conditions 


(o  D“-*jt)(0+)  =  c*.  k=  1 . n.  (19) 

The  notation  (o D“_*jc)(0+)  means  that  the  limit  is  taken  at  almost  all  points  of  the  right-sided 
neighborhood  (0,  0  +  e),  e  >  0,  of  zero  as  follows: 


(0D“~kx)(0+)  =  Km  oD“~kx(t),  (k  =  1 . n  -  1). 

(0D“-"jt)(0+)  = 


The  Cauchy-type  problem  (18)  and  (19)  can  be  reduced  to  the  nonlinear  Volterra  integral 
equation  of  second  kind 


f«)  =  ]C 


Ck 


k=\ 


r  (a  —  k  +  1) 


1  r1  G[r,jt(r)]dr 
F(a)  Jo  (t  -  r)1_a 


(20) 


where  t  >  0.  The  result  was  obtained  by  Kilbas,  Bonilla  and  Trujillo  in  [21,  22].  For 
a  =  n  =  2,  equation  (20)  gives  (10). 

The  Cauchy-type  problem  (18)  and  (19)  and  the  nonlinear  Volterra  integral  equation  (20) 
are  equivalent  in  the  sense  that,  if  x(t)  e  L(  0,  tf)  satisfies  one  of  these  relations,  then  it  also 
satisfies  the  other.  In  [21,  22]  (see  also  theorem  3.1.  in  section  3.2  1  of  [7]),  this  result  is 
proved  by  assuming  that  the  function  G[f,  x]  belongs  to  L( 0,  t/)  for  any  x  e  W  c  R. 

Let  us  give  the  basic  theorem  regarding  the  nonlinear  differential  equation  involving  the 
Riemann-Liouville  fractional  derivative. 


Kilbas-Bonilla-Trujillo  theorem.  Let  W  be  an  open  set  in  R  and  let  G  [f ,  x  ],  where  t  €  (0,  //  ] 
and  x  €  W,  be  a  real-valued  function  such  that  G[r,  x]  €  L(0,  tf)  for  any  x  €  W.  Let  x(t) 
be  a  Lebesgue  measurable  function  on  (0,  t /).  If  x(t)  6  L(  0,  tf),  then  x(t)  satisfies  almost 
everywhere  equation  (18)  and  conditions  (19)  if  and  only  if  x(t)  satisfies  almost  everywhere 
the  integral  equation  (20). 


Proof.  This  theorem  is  proved  in  [21,  22]  (see  also  theorem  3.1.  in  section  3.2.1  of  [7]).  □ 


In  [14]  we  consider  a  fractional  generalization  of  equation  (1)  of  the  form 

0D^(0  +  *G[A:(f)]fS(^-*)  =0.  (1  <  a  <  2),  (21) 

k=\  '  ' 

where  t  >  0,  and  o D*  is  the  Riemann-Liouville  fractional  derivative  defined  by  (17).  Let  us 
give  the  following  theorem  for  equation  (21). 
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Proposition  2.  The  Cauchy-type  problem  for  the  fractional  differential  equation  of  the  form 
(21)  with  the  initial  conditions 

(0D“-1ac)(0+)=c1,  (gD“~2x)(0+)  =  (0/,2-“x)(0+)  =  c2  (22) 

is  equivalent  to  the  equation 


x(t)  -  ta  1  +  • 


ci 


r(a  -  1) 

where  nT  <  t  <  (n  +  \  )T. 


ta~ 2  - 


KT 


r(“)£r 


£  G[*(*7~)](/-*7y-\ 


(23) 


Proof.  Using  the  function 

G[t,x(0]  =  -KG[x]Yjs(j-kY  (24) 

*=1  '  ' 

equation  (21)  has  the  form  of  (18)  with  the  Riemann-Liouville  fractional  derivative  of  order 
a,  where  1  <  a.  ^  2.  It  allows  us  to  use  the  Kilbas-Bonilla-Trujillo  theorem.  As  a  result, 
equation  (2 1 )  with  initial  conditions  ( 1 9)  of  the  form  (22)  is  equivalent  to  the  nonlinear  Volterra 
integral  equation 


*(/)  = 


r» 


ria^Tj'0"2  -  rS  gi  d'  c s  (f  - l)  • 


(25) 


where  /  >  0.  If  nT  <  t  <  (n  +  1)7,  then  the  integration  in  (25)  with  respect  to  r  gives  (23). 
This  ends  the  proof.  □ 


To  obtain  equations  of  discrete  map  a  momentum  must  be  defined.  There  are  two 
possibilities  of  defining  the  momentum: 

p(0  =  0  Dp'xU),  p(t)  =  D)x(t).  (26) 

Let  us  use  the  first  definition.  Then  the  momentum  is  defined  by  the  fractional  derivative 
of  order  a  —  1 .  Using  the  definition  of  the  Riemann-Liouville  fractional  derivative  (17)  in  the 
form 

0D?x(t)  =  D1,  oI?~ax(t),  (1<«sS2),  (27) 

we  define  the  momentum 

1  Id  f*  ;t(r)dr 

P(t)  =  oD“~lx(t)  =  p^Ta)  d/  Jo  (T-tF"*  ’  (,<^2)’  (28) 

where  *(r)  is  defined  for  r  e  (0,  /).  Then 

0  D?x(t)  =  Djp(t),  (1  <a^2).  (29) 

Using  momentum  p(t)  and  coordinate  x(t ),  equation  (21)  can  be  represented  in  the 
Hamiltonian  form 


0  D?-'x(t)  =  p(t),  (30) 

D}p(t)  =  -KG[x(t)]jfi&U;-ky  (31) 
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Proposition  3.  The  Cauchy-type  problem  for  the  fractional  differential  equations  of  the  form 
(30)  and  (31)  with  the  initial  conditions 

(0D“-'*)(0+)=Ci,  (o/J>r2*)(0+)  =  (o/,2-"*)(0+)  =  Q  (32) 

is  equivalent  to  the  discrete  map  equations 


xn+\  = 


c\Ta 


-1 


r» 


-(n  +  1)"  1  + 


r^Ta~^  KTa  _n 

:0l  +  l)<‘-2___y2C[xk}(n  +  \-kr 


!>-!) 


-I 


r(a) 


(33) 


pn+ 1  =  ci  -  at  yy  c[x*]. 

k=i 


(34) 


Proof.  We  use  proposition  2  to  prove  this  statement.  To  obtain  an  equation  for  the  momentum 
(28),  we  use  the  following  fractional  derivatives  of  power  functions  (see  section  2.1  in  [7]): 

aD f(r  -  a/-1  =  — ■■  (?  -  a)*-1-".  a  S*  0,  0  >  0,  t  >  a,  (35) 

l(p  —  a) 

0D“r“_*=0,  *  =  1 . n,  « —  1  <a^n.  (36) 

These  equations  give 

oDf*"-1  =  r(a).  0O(“r“-2  =  0 

and 

aD«{t-a)°-x  =  r(or). 

We  note  that  equation  (23)  for  x(t)  can  be  used  only  if  t  €  ( nT ,  /),  where  /i7  <  t  < 
(n  +  1)7.  The  function  ^(r)  in  the  fractional  derivative  oD“  of  the  form  (28)  must  be  defined 
for  all  r  €  (0,  /).  We  cannot  take  the  derivative  0 Df  of  the  functions  (r  -  kT)a~l  that  are 
defined  for  r  e  ( kT ,  t).  In  order  to  use  equation  (23)  on  the  interval  (0,  r),  we  must  modify 
the  sum  in  equation  (23)  by  using  the  Heaviside  step  function.  Then  equation  (23)  has  the 
form 

* (t)  =  tt"-'  +  -  ra~2  -ffrY  G[x(kT)](T-kT)a-'9(T-kT),  (37) 

r(a)  r(a-i)  r(o;)  " 

where  r  e  (0,  r).  Using  the  relation 

0D“ ((f  —  a)“~l  0{t  —  a))  —  aD?(t -a)a~l  =  r(a),  (38) 

equations  (28)  and  (37)  give 

n 

Pit)  =  c,  -  KT  J2  G[x(.kT)],  (39) 

k=  1 

where  nT  <  t  <  (n  +  1)7.  Then  the  solution  of  the  left  side  of  the  (n  +  1 ) th  kick 

n 

Pn+t  =  c,  -  KT  G[JC*3.  (40) 

*~1 

As  a  result,  we  obtain  a  universal  map  with  memory  in  the  form  of  equations  (33) 
and  (34). 

This  ends  the  proof.  □ 
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Remark  3.  For  a  =  n  =  2  equations  (33)  and  (34)  give  the  usual  universal  map  (6)  and  (7). 
Remark  4.  We  note  that  the  map  (33)  and  (34)  with 
C\  =  P\,  C2  =  0 
was  obtained  in  [14]  in  the  form 

-*/i+i  =  L—  53  PMVa(n  - k  +  1),  (41) 

P,*i  =Pn~  KTG(xn),  (1  <  a  <  2),  (42) 


where  p\  =  c\ ,  and  the  function  Va(z)  is  defined  by 

Va(z)  =  z“_1  —  (z  —  1)“-1,  (z^l).  (43) 

In  [14],  we  obtain  these  map  equations  by  using  an  auxiliary  variable  f  (f)  such  that 

£D,2-“£(0  =  x(r). 

The  nonlinear  Volterra  integral  equations  and  the  general  initial  conditions  (32)  are  not  used 
in  [14].  In  the  general  case,  the  fractional  differential  equation  of  the  kicked  system  (21)  is 
equivalent  to  the  discrete  map  equations 


Xn+1  =  7T— tT!  PMVa(n  -  k  +  \)  +  — —  (n  +  l)“-2,  (44) 

r(ar)  “  T(a-l) 

p„+1  =pn-  KTG(xn ),  (1  <  a  2),  (45) 


where  p\  =  C\.  Here  we  take  into  account  the  initial  conditions  (32).  The  second  term  of  the 
right-hand  side  of  equation  (44)  is  not  considered  in  [14].  Using  —1  <  a  —  2  <  0,  we  have 

lim  (n  +  1)“-2  =  0. 

rt— ►OO 

Therefore,  the  case  of  large  values  of  n  is  equivalent  to  C2  =  0. 

Let  us  give  the  proposition  regarding  the  second  definition  of  the  momentum  p(t)  = 
Djx(t). 

Proposition  4.  The  Cauchy-type  problem  for  the  fractional  differential  equations 


Djx(t)  =  p(t). 

(46) 

0  D?x(t)  =  -KG[x(t)]  J2  8  (j;  -  ^  , 

(1  <  a  <  2) 

(47) 

with  the  initial  conditions 

(oD“-^)(0+)  =  ci,  (qD°-2x)(0+)  =  (0/2-“at)(0+)  =  c2  (48) 

is  equivalent  to  the  discrete  map  equations 


C]  7’“-' 

r(a) 


(«  +  !)“-' 


+ 


c2r°-2 
T(a  -  1) 


(n  +  1)“-2 


KTa 

T>) 


y~l  GUt](n  +  1  -k)a  \ 

k=  1 


(49) 
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c  7-0-2 
P„+,  =  =7 - 7T('I  +  »“ 


r(a-  1) 


2  C2(a  —  2)Ta  3  , 

+  - — — (n  +  1)“~3 


r(a  -  1) 


Kja-l  n 


a -2 


(50) 


Proof.  We  define  the  momentum 
pit)  =  Djx(t). 

If  nT  <  t  <  (n  +  1)7,  then  the  differentiation  of  (23)  with  respect  to  t  gives 


pit)  =  +  £il£L^,«-3  _  _SL_ 

1  r(a-l)  r(a-l)  r(a-l) 


—  G[x(kT)](t  -  kT)a~2. 


(51) 


k- 1 


Here  we  use  the  relation 

r(a)  =  (a  -  l)T(a  -  1),  (1  <  a  ^  2). 

Using  equations  (23)  and  (51),  we  can  obtain  the  solution  of  the  left  side  of  the  (a  +  l)th  kick 
(13)  and  (14).  As  a  result,  we  have  equations  (49)  and  (50). 

This  ends  the  proof.  □ 


Remark  5.  Equations  (49)  and  (50)  describe  a  generalization  of  equations  (6)  and  (7).  If 
a  =  n  =  2  and  ci  =  jcq,  c\  =  po»  then  equation  (49)  gives  (6)  and  (7). 


Remark  6.  In  equations  (50)  and  (5 1 ),  we  can  use 

C2(a  ~  2)  _  c2 
T(a  -  1)  ”  T(a  — 2) 

for  1  <  a  <  2. 


Remark  7.  If  we  use  the  definition  p(t)  =  Djx(t)y  then  the  Hamiltonian  form  of  the 
equations  of  motion  will  be  more  complicated  than  (30)  and  (31)  since 

D2Ql}-ax{t)  #  or--aD;x(t). 


Remark  8.  Note  that  we  use  the  usual  momentum  p(t)  =  Djx(t).  In  this  case,  the  values 
cj  and  C'i  arc  not  connected  with  p( 0)  and  jt(0).  If  we  use  the  momentum  p(t)  =  oD^lx(t)y 
then  ci  =  p( 0). 


4.  Riemann-Liouville  and  Caputo  fractional  derivatives 

In  [14]  wc  consider  nonlinear  differential  equations  with  Riemann-Liouville  fractional 
derivatives.  The  discrete  maps  with  memory  are  obtained  from  these  equations.  The  problems 
with  initial  conditions  for  the  Riemann-Liouville  fractional  derivative  are  not  discussed. 

The  Riemann-Liouville  fractional  derivative  has  some  notable  disadvantages  in 
applications  in  mechanics  such  as  the  hyper-singular  improper  integral,  where  the  order  of 
singularity  is  higher  than  the  dimension,  and  nonzero  of  the  fractional  derivative  of  constants, 
which  would  entail  that  dissipation  docs  not  vanish  for  a  system  in  equilibrium.  The  desire 
to  use  the  usual  initial  value  problems  for  mechanical  systems  leads  to  the  use  of  Caputo 
fractional  derivatives  [7,  6]  rather  than  the  Riemann-Liouville  fractional  derivative. 
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The  left-sided  Caputo  fractional  derivative  [7,  23-25]  of  order  a  >  0  is  defined  by 


oC0,"/(0  = 


0 /,"/(')  = 


r(n-a)  JQ 


dr  D;/( r) 

(f  -  r)a~n^ 

where  n  —  \  <  a  <  n,  and  o /“  is  the  left-sided  Riemann-Liouville  fractional  integral  of  order 
a  >  0  that  is  defined  by 

This  definition  is,  of  course,  more  restrictive  than  the  Riemann-Liouville  fractional 
derivative  [4,  7]  in  that  it  requires  the  absolute  integrability  of  the  derivative  of  order  n. 
The  Caputo  fractional  derivative  first  computes  an  ordinary  derivative  followed  by  a  fractional 
integral  to  achieve  the  desire  order  of  fractional  derivative.  The  Riemann-Liouville  fractional 
derivative  is  computed  in  the  reverse  order.  Integration  by  part  of  (52)  will  lead  to 


r  /( r)dt 
Jo  (/  —  r)1-"’ 


C t  >  0). 


(52) 


(53) 


c0D°x{t)  =  oD7*<0  -  E  r7r-77nx<t)ffl>-  (54) 

It  is  observed  that  the  second  term  in  equation  (54)  regularizes  the  Caputo  fractional  derivative 
to  avoid  the  potentially  divergence  from  singular  integration  at  /  =  0.  In  addition,  the  Caputo 
fractional  differentiation  of  a  constant  results  in  zero 

o  D“C  =  0. 

Note  that  the  Riemann-Liouville  fractional  derivative  of  a  constant  need  not  be  zero  and  we 
have 


o 


D?C  = 


ra 

r(i-a) 


c. 


If  the  Caputo  fractional  derivative  is  used  instead  of  the  Riemann-Liouville  fractional 
derivative,  then  the  initial  conditions  for  fractional  dynamical  systems  are  the  same  as  those 
for  the  usual  dynamical  systems.  The  Caputo  formulation  of  fractional  calculus  can  be  more 
applicable  in  mechanics  than  the  Riemann-Liouville  formulation. 


5.  Caputo  fractional  derivative  and  universal  map  with  memory 

In  this  section,  we  study  a  generalization  of  differential  equation  (1)  by  the  Caputo  fractional 
derivative.  The  universal  map  with  memory  is  derived  from  this  fractional  equation 

We  consider  the  nonlinear  differential  equation  of  order  or,  where  0  ^  n  —  1  <  a  ^  n, 

?0*x(f)-G[r.x(O).  (0  <»<</).  (55) 

involving  the  Caputo  fractional  derivative  D“  on  a  finite  interval  [0,  t /  ]  of  the  real  axi  s,  with 
the  initial  conditions 

(D*x)(0)  =  c*,  k  =  0, . . . ,  n  —  1 .  (56) 

Kilbas  and  Marzan  [19,  20]  proved  the  equivalence  of  the  Cauchy-type  problem  of  the  form 
(55),  (56)  and  the  Volterra  integral  equation  of  second  kind 

*(0  =  E  tS  +  F7T  /  dr  Gtr-  *<*>]  (t  ~  r)a~'  (57) 

tok'  r(a) 

in  the  space  Cn~l  [0,  t/].  For  a  =  n  =  2  equation  (57)  gives  (10). 
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Let  us  give  the  basic  theorem  regarding  the  nonlinear  differential  equation  involving  the 
Caputo  fractional  derivative. 

Kilbas-Marzan  theorem.  The  Cauchy-type  problem  (55)  and  (56)  and  the  nonlinear  Volterra 
integral  equation  (57)  are  equivalent  in  the  sense  that,  if  x(t)  €  C[0,  tf]  satisfies  one  of  these 
relations ;  then  it  also  satisfies  the  other. 

Proof.  In  [19,  20]  (see  also  [7],  theorem  3.24.)  this  theorem  is  proved  by  assuming  that  a 
function  G[/,jc]  for  any  jc  6  W  C  R  belongs  to  CK(0, //)  with  0  ^  y  <  1,  y  <  ct.  Here 
C/(0,  tf)  is  the  weighted  space  of  functions  /[/]  given  on  (0,  //],  such  that  tY  f[t]  €  C(0,  tf). 


This  ends  the  proof.  □ 

We  consider  the  fractional  differential  equation  of  the  form 

<-D?xO)  +  KG[x(t)]J2s(  Y~k)  =0-  (1  <  a  <  2),  (58) 

k=\  '  ' 

where  £  Da{  is  the  Caputo  fractional  derivative,  with  the  initial  conditions 

Jf(0)  =  JCo,  (D'jc)(0)  =  po  (59) 

Using  p{t)  =  D]x(t)%  equation  (58)  can  be  rewritten  in  the  Hamilton  form. 

Proposition  5.  The  Cauchy-type  problem  for  the  fractional  differential  equations 

D}x(t)  =  p(t),  (60) 

£z>rV(«)  =  -tfG[*(»)]£5(f -*J.  (1  <  a  <  2)  (61) 


k- 1  X  7 

with  the  initial  conditions 

x(Q)  =  x0,  p(0)  =  p0 

is  equivalent  to  the  discrete  map  equations 

A:n+i  =x0  +  p0(n  +  1)7  -  — —  Y  (n  +  1  -*)"“' C[ac*], 

Pn  +  i  =  Po-  -pT - 77  Y  («  +  •  -  k)a~2G[.Xk]. 

r(ar  -  I) 

Proof.  Wc  use  the  Kilbas-Marzan  theorem  with  the  function 
G[;,*(0]  = -KG[x(;)J  ^ 

k  =  1  '  7 

The  Cauchy-type  problem  (58)  and  (59)  is  equivalent  to  the  Volterra  integral  equation  of 
second  kind 

x(O  =  x0+p0t-d^'y2j  dr  (t  -  r)a~l  G[x(r)]S  -  k^j  ,  (65) 

in  the  space  of  continuously  differentiable  functions  x(t)  €  C'[0,  tf] 


(62) 

(63) 

(64) 
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If  nT  <  t  <  (n  +  \)T,  then  equation  (65)  gives 


(66) 


p(t)  =  Djx(t). 

Then  equations  (66)  and  (67)  give 


(67) 


(68) 


where  we  use  T (a)  =  (a  —  1 ) T  (a  —  1 ). 

The  solution  of  the  left  side  of  the  ( n  +  l)th  kick  (13)  and  (14)  can  be  represented  by 
equations  (63)  and  (64),  w  here  we  use  the  condition  of  continuity  x(tn  +  0)  =  x(tn  —  0). 


This  ends  the  proof. 


□ 


Remark  9.  Equations  (63)  and  (64)  define  a  generalization  of  the  universal  map.  This 
map  is  derived  from  a  fractional  differential  equation  with  Caputo  derivatives  without  any 
approximations.  The  main  property  of  the  suggested  map  is  a  long-term  memory  that  means 
that  their  present  state  depends  on  all  past  states  with  a  power-law  form  of  weights. 

Remark  10.  If  or  =  2,  then  equations  (63)  and  (64)  give  the  universal  map  of  the  form  (6) 
and  (7)  that  is  equivalent  to  equations  (2).  As  a  result,  the  usual  universal  map  is  a  special 
case  of  this  universal  map  with  memory. 

Remark  11.  By  analogy  with  proposition  5,  it  is  easy  to  obtain  the  universal  map  with 
memory  from  fractional  equation  (58)  with  or  >  2. 

6.  Conclusion 

The  suggested  discrete  maps  with  memory  are  generalizations  of  the  universal  map.  These 
maps  describe  fractional  dynamics  of  complex  physical  systems.  The  suggested  universal 
maps  with  memory  are  equivalent  to  the  correspondent  fractional  kicked  differential  equations. 
We  obtain  a  discrete  map  from  a  fractional  differential  equation  by  using  the  equivalence  of 
the  Cauchy-type  problem  and  the  nonlinear  Volterra  integral  equation  of  second  kind.  An 
approximation  for  fractional  derivatives  of  these  equations  is  not  used. 

It  is  important  to  obtain  and  to  study  discrete  maps  which  correspond  to  the  real  physical 
systems  described  by  the  fractional  differential  equations.  In  mechanics  and  electrodynamics, 
we  can  consider  viscoelastic  and  dielectric  materials  as  media  with  memory.  We  note  that  the 
dielectric  susceptibility  of  a  wide  class  of  dielectric  materials  follows,  over  extended  frequency 
ranges,  a  fractional  power-law  frequency  dependence  that  is  called  the  ‘universal’  response 
[15,  16].  As  was  proved  in  [17,  18],  the  electromagnetic  fields  in  such  dielectric  media 
are  described  by  differential  equations  with  fractional  time  derivatives.  These  fractional 
equations  for  electromagnetic  waves  in  dielectric  media  are  common  to  a  wide  class  of 
materials,  regardless  of  the  type  of  physical  structure,  chemical  composition,  or  of  the  nature 
of  the  polarizing  species,  whether  dipoles,  electrons  or  ions.  We  hope  that  it  is  possible  to 
obtain  the  discrete  maps  with  memory  which  correspond  to  the  real  dielectric  media  described 
by  the  fractional  differential  equations. 
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Abstract 

Using  kicked  differential  equations  of  motion  with  derivatives  of  non-integer  orders,  we  obtain  gen¬ 
eralizations  of  the  dissipative  standard  map.  The  main  property  of  these  generalized  maps,  which  are 
called  fractional  maps,  is  long-term  memory.  The  memory  effect  in  the  fractional  maps  means  that 
their  present  state  of  evolution  depends  on  all  past  states  with  special  forms  of  weights.  The  fractional 
dissipative  standard  maps  are  used  to  study  attractors  of  the  systems  described  by  kicked  fractional 
differential  equations. 
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Discrete  maps  are  widely  used  to  study  general  properties  of  dynamical  sys¬ 
tems.  In  those  cases  when  they  can  be  derived  from  differential  equations,  their 
analysis  gives  the  exact  properties  of  the  corresponding  systems.  In  this  ar¬ 
ticle  we  derive  discrete  maps  (fractional  maps)  from  the  fractional  differential 
equations,  which  correspond  to  the  fractional  generalizations  of  the  dissipative 
standard  map  [1].  We  demonstrate  how  the  attractors  of  the  fractional  maps 
are  different  from  the  attractors  of  the  dissipative  standard  map. 

I.  INTRODUCTION 

There  is  a  number  of  distinct  areas  of  physics  where  basic  problems  can  be  reduced  to 
the  study  of  simple  discrete  maps.  Discrete  maps  as  substitutes  of  differential  equations 
have  been  used  to  study  evolution  problems  in  [2-6].  They  lead  to  a  simpler  formalism, 
which  is  particularly  useful  in  simulations.  The  dissipative  standard  map  is  one  of  the  most 
widely  studied  maps.  In  this  paper  we  consider  fractional  generalizations  of  the  dissipative 
standard  map  which  are  described  by  fractional  differential  equations  [7  9]. 

The  treatment  of  nonlinear  dynamics  in  terms  of  discrete  maps  is  a  very  important  step  in 
understanding  the  qualitative  behavior  of  systems  described  by  differential  equations.  The 
derivatives  of  non-integer  orders  are  a  natural  generalization  of  the  ordinary  differentiation 
of  the  integer  order.  The  fractional  differentiation  with  respect  to  time  is  characterized 
by  long-term  memory  effects  which  correspond  to  intrinsic  dissipative  processes  in  physical 
systems.  The  application  of  memory  effects  to  discrete  maps  means  that  their  present  state 
of  evolution  depends  on  all  past  states  [10-15,  17]. 

Discrete  maps  with  memory  can  be  derived  (see  [17])  from  equations  of  motion  with 
fractional  derivatives.  In  Ref.  (17  a  fractional  generalization  of  the  standard  map  has  been 
derived  from  a  fractional  differential  equation.  A  fractional  generalization  of  the  dissipative 
standard  map  was  also  suggested  in  [17  .  Unfortunately,  in  that  generalization  a  dissipation 
was  introduced  by  the  change  of  the  variable  pn  ->  — bpn .  The  map  equations  were  not 
directly  connected  with  a  fractional  equation  of  motion.  In  this  paper  we  propose  two 
generalizations  of  the  dissipative  standard  map.  The  first  one  is  derived  from  a  differential 
equation  with  fractional  damped  kicks.  The  second  generalization  of  the  dissipative  standard 
map  is  derived  from  a  fractional  differential  equation  (kickes  are  not  fractional)  .  A  nonlinear 
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system  with  fractional  derivatives  perturbed  by  a  periodic  force  exhibits  a  new  type  of  chaotic 
motion  which  can  be  called  the  fractional  chaotic  attractor.  Fractional  discrete  maps  [17]  are 
used  to  study  new  types  of  attractors  of  fractional  dynamics  described  by  kicked  fractional 
equations.  In  this  paper  some  fractional  differential  equations  of  motion  of  kicked  systems 
with  friction  are  considered.  Corresponding  discrete  maps  with  memory  are  derived  from 
these  equations.  The  fractional  generalizations  of  the  dissipative  standard  map  are  suggested 
and  these  maps  are  used  in  computer  simulations. 


II.  DISCRETE  MAPS  WITHOUT  MEMORY 

In  this  section,  a  brief  review  of  discrete  maps  is  considered  to  fix  notations  and  provide 
convenient  references.  For  details,  see  [2-6] . 


A.  Standard  map 


Let  us  consider  the  equation  of  motion 


oc 

K  siii(x)  ^  s(t  —  nj  =  0, 


n~0 


(1) 


in  which  perturbation  is  a  periodic  sequence  of  delta-function  type 
with  period  T  =  1.  A'  is  the  amplitude  of  the  pulses.  This  equation 
Hamiltonian  form 


OC 

i  =  p,  |i  + A'  sin(x)  -  rij  =  0. 


n~0 


pulses  (kicks)  following 
can  be  presented  in  the 

(2) 


It  is  well-known  that  these  equations  can  be  represented  (see  for  example  Chapter  5  in  [3]) 
in  the  form  of  discrete  map 


•Tn-H  -Tn  "1“  Pn+1* 


Pn41  =  Pn  ~  A  Sin(Xn). 


(3) 

(4) 


Equations  (3)  and  (4)  are  called  the  standard  map.  This  map  is  also  called  the  Chirikov 
map  (4j. 
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B.  Dissipative  standard  map 


The  dissipative  standard  rnap  [1,  18,  19]  is 


1 1  =  Xn  +  (iYn+ i  +  a 

(5) 

K.  +  1  =  ('~q  [>n  +  e  sin(A'n)], 

(6) 

where  jj,  =  (eq  —  1)/^,  The  dissipative  standard  map  is  also  called  the  Zaslavsky  rnap.  Note 
that  a  shift  Q  doevS  not  play  an  important  role  and  it  can  be  put  to  aero  (Q  =  0).  The 
dissipative  standard  map  with  Q  =  0  can  be  represented  by  the  equations 

-Xvii  =  Xn  +  Pn+ (7) 

Pn+ 1  =  ~bPn  -  Z  fAn(Xn).  (8) 

For  the  parameters 

Z  =  Pn  =  nYn,  6=-fT*  (9) 

equations  (7)  and  (8)  give  Eqs.  (5)  and  (6)  with  il  —  0. 

For  b  =  —  1  and  Z  =  A\  we  get  the  standard  map  which  is  described  by  Eqs.  (3)  and  (4) 
with  T  =  1 . 

Note  that  for  large  q  -+  oc  (for  small  h  -»  0)  Eqs.  (7)  and  (8)  with  Z  =  ~A'  shrink  to 
the  proposed  by  Arnold  [20  one-dimensional  sine-map 

Xn+i**Xn  +  K*m(Xn).  (10) 


C.  Kicked  damped  rotator  map 

The  equation  of  motion  for  kicked  damped  rotator  is 

x  -r  qx  —  KG(x)  ^  S  —  nT ^ . 

tl'ssK) 

It  is  well-known  [5]  that  Eq.  (11)  gives  the  two-dimensional  map 

1  —  c  ^ 

^n+l  :=  *  d  “  [Pn  “F  JKG{xri  )^ , 

Pn  + 1  =  C",r[Pn  +  ^G(*n);. 


(11) 


(12) 

(13) 
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This  map  is  known  as  the  kicked  damped  rotator  rnap*  The  phase  volume  shrinks  each  time 
step  by  a  factor  cxp(— 7).  The  map  is  defined  by  two  important  parameters,  dissipation 
constant  7  and  force  amplitude  K.  These  equations  can  be  rewritten  in  the  form 

eqI  -  1 

n  }  i  3'n  T  “  Pni  1 i 
Pn,  \  =  e~ql  [p„  +  KG{xn)}. 

It  is  easy  to  see  that  these  equations  give  the  dissipative  standard  map  (7)  mid  (8)  with 
fi  =  0  if  we  use 

X„  =  xn.  Yn  =  pn.  e  =  A',  T=  1,  G'(j-)  =  sin(j:) 

This  allows  us  to  derive  dissipative  standard  map  (5)  and  (6)  from  the  differential  equation. 

X  +qX  =eain(X)J2s(t-n)-  (14) 

n~  0 

These  equations  give  the  discrete  map  defined  by  Eqs.  (5)  and  (6)  with  Q  =  0. 


III.  FRACTIONAL  STANDARD  MAP  AND  DISSIPATION 


A  fractional  generalization  of  the  differential  equation  (1)  has  been  suggested  in  [17^.  The 
discrete  map  which  corresponds  to  the  fractional  equation  of  order  1  <  o  <  2  was  derived. 
This  map  can  be  considered  as  a  generalizat  ion  of  the  standard  map  for  the  case  1  <  a  <  2. 
We  consider  a  fractional  generalization  of  (1)  in  the  form 


ac 

o Dfx  +  K sin(r)  ^ ~  =  0.  (1  <  a  <  2), 


n=0 


where  o D*  is  the  Riemann-Lioiiville  fractional  derivative  [7-9],  which  is  defined  by 
n»T  _  n 2  j2  a  1  ^  f  :E(T)^'  fl  <r  o  <  21 

0Dtx  Dt  0/,.  -r  p(2  —  a)  dP  JQ  (t  —  t)“-1  ’ 

Here  we  use  the  notation  D2,  =  rP /dt2,  and  0/“  is  a  fractional  integration  [7-9]. 
Defining  the  momentum  as 

p(t)  =  o  D?~lx(t), 


(15) 


(16) 


and  using  the  initial  conditions 


(OJDr1n:)(0+)  =  Pu  (oDr2x)(0+)  -  b , 


(17) 


o 


it  is  possible  to  derive  the  equation  for  the  fractional  standard  map. 


Proposition  1.  The  fractional  differential  equation  of  the  kicked  system  (15)  is  equiva¬ 
lent  to  the  discrete  map 

-Tn+1  =  r(£)  ~  k  + !)  +  +  1)a'2?  (18) 

Pw+1  =  Pn  -  K  sin(xn),  (1  <  a  <  2 ),  (19) 

where  the  function  Va(z)  is  defined  by 

Va(z)  =  ^  -  (z  -  1)" . l.  (20) 

Proof  of  this  Proposition  is  given  in  [21 

A  fractional  generalization  of  the  dissipative  standard  map  suggested  in  [1,  18  can  be 
defined  by 

1  n 

*«+i  =  Yj-)22pkMn-k  +  \),  (1  <  a  <  2),  (21) 

Pnei  =  -frPn  -  z  sm{xn)<  (22) 

where  the  parameters  are  defined  by  conditions  (9).  For  h  =  —  1  and  Z  =  K  Eqs.  (21) 
and  (22)  give  the  fractional  standard  map  with  T  =  1.  Note  that  this  fractional  dissipative 
standard  map  is  not  derived  from  a  fractional  differential  equation.  This  map  is  derived 
by  pn  —  fepn  in  the  fractional  standard  map.  Fractional  dissipative  standard  map  can  be 
derived  from  fractional  differential  equations.  In  this  paper,  we  derive  two  fractional  gen¬ 
eralizations  of  the  dissipative  standard  map  which  are  obtained  from  fractional  differential 
equations. 

IV.  FRACTIONAL  DERIVATIVE  IN  THE  KICKED  TERM  AND  THE  FIRST 
FRACTIONAL  DISSIPA1IVE  STANDARD  MAP 

In  this  section  we  suggest  the  first;  fractional  generalization  of  differential  equation  (11) 
for  a  kicked  damped  rotator.  In  this  generalization  we  introduce  a  fractional  derivative  in 
the  kicked  damped  term,  i.e.  the  term  of  a  periodic  sequence  of  delta-function  type  pulses 
(kicks),  and  derive  the  corresponding  discrete  map. 


6 


Consider  the  fractional  generalization  of  equation  (11)  in  the  form 


oc 

D2t  X{l)  -  qD\X{t)  =  c  sin  (£Z>fx)  JZ  6(l  ~ Ti)>  (°  </?<!), 


(23) 


n— 0 


where  g  €  R,  and  „  Df  is  the  Caputo  fractional  derivative  [9;  of  the  order  0  <  B  <  l  defined 
bv 

(24) 

Here  we  use  the  notation  D\  —  d/dt)  and  q /?  is  a  fractional  integration  [7-9; .  For  0  =  0 
fractional  equation  (23)  gives  equation  (11).  Note  that  we  use  the  minus  on  the  left-hand 
side  of  Eq.  (23),  where  q  can  be  a  positive  or  negative  value.  Fractional  derivative  q  DtX 
is  presented  in  the  kicked  damped  term. 


Proposition  2,  The  fractional  differential  equation  of  the  kicked  system  (23)  is  equiva¬ 
lent  to  the  dis  erode  map 

i  __  «-<? 

(25) 


XnH  i  -  +  — — W 


<7 

n~  1 


v*+ 1  =  yn  +  £  sin(f[Y~y  5Z  Yk+t  k  - «)) , 

where  the  functions  b)  are  defined  by 

Wa_fl(a,h)  =  a'1  lea(b+r>  F(1  -  ft, ah)  -  P(1  -  B,a{h  +  1)) 
and  r(n,  b)  is  the  incomplete  Gamma  function 

r(o,  h)  —  J(°  ya~le~vdy. 

Proof.  Fractional  equation  (23)  can  he  presented  in  the  Hamiltonian  form 

X  =  Y , 

30 

Y  -  qY  =  c  sin(  $D?X)  ^  <*(*  ~  n)- 


n~Q 


where  0  <  3  <  0,  and  q  €  R. 
Between  any  two  kicks 


Y  -  qY  «  0. 


(26) 


(27) 


(28) 


7 


(29) 


For  L  £  ( t.n  +  0,  l.n ( i  —  0),  the  solution  of  Eq.  (29)  is 

V'ftn+l  -0)  =  Y(tn+Q)C<l.  (30) 

Let  us  use  the  notations  tn  =  nT,  with  T  =  1  and 


A'„  =  X(tn  —  0)  =  limX(n  —  r), 

e~H) 


Yn  =  Y(tn  —  0)  =  lim  r(n  —  f) 

f  —►(} 


For  t  €  (in  —  (Jn+i  —f),  the  general  solution  of  (28)  is 

Y(t)  =  Yneq[t  £"}  +  e  ]T  M o  D?mX)  f  dTt/*t-T)6(T  -  m). 


Then 


rn+1  —  ea  jVn  +  e  sin(  c0DfnX) 
Using  (33),  the  integration  of  the  first  equation  of  (28)  gives 


*»  +  !-*» 


1  -e« 


Yn+eMoDtX) 


(31) 


(32) 


(33) 


(34) 


Let  us  consider  the  Caputo  fractional  derivative  from  Eqs.  (33)  arid  (34).  It  is  defined  by 
the  equation 


=0 


l-.3n  l 


rtrmf. 


dr  dX{r) 


l'(l  -  .3)  J0  \U  -Tf  dr 

Using  Y (r)  =  dX  {r)/dr.  this  relation  can  be  rewritten  as 


,  (O<0<1). 


C  nJ  y  _  1  V'  ["»  Y(T)dT 

0  £"  r(i  u,  -  ry' 


(35) 


where  £fc+i  =  t.k  +  1  =  (k  +  1).  and  tk  —  k,  such  that  to  =  0.  For  r  £  (tk.tk+i),  equations 
(30)  and  (31)  give 

Y{ t)  =  Y(tk  +  0)e?(r  =  F(4+1  -  0)e  = 

=  Yk^e'Y  “-l)  =  Yk+le<Y-^. 


Then 


Jdi 
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Yk+l  f  "  *  t* *•-*>»— )z-*dz  =  n+1e?(t"-£‘*l)  f"  z~^e~q*dz 


rqitn  ~  tfc ) 

=  y*+1/  /  y~3e~vdy. 

JqU 


As  a  result,  equation  (36)  gives 


/t/(en  -e^j) 

rtk"  Y(r)dr 


(36) 


rtk"  r(r 

Jtk  {tn-T)~?- 

fr(i  -  A  <?(*„  -  4+i))  -  r(i  -  0,  q(tn  -  tfc)) 


(37) 


Here  F(a,  6)  is  the  incomplete  Gamma  function  (27),  where  a  and  h  are  complex  numbers. 
Using  (35)  and  (37),  we  obtain 

n-l 


oKx 


i 


F(1  -  8) 


Y,Y^W2„3(<hk~n).  (0  <  3  <  1), 


k=t 0 


where 


W9-p{at  b)  =  l)  F(1  -  A  ab)  -  F(1  -  ,5,  a{b  +  1)) 

Substitution  of  (38}  into  (33)  and  (3-4)  gives 

n  } 

Yn+ 1  =  c* [y„  +  esin(--  i . —  £  WF2 -*(<?•  *  -  n)) 


Vn+i  —  Xn - - —  Ki  +  £  sin(— r  ^  y*+xW/2-^(?>  ~  n)) 

Equations  (40)  and  (41)  can  be  presented  in  the  form  of  Eqs.  (25)  and  (26). 
This  ends  the  proof.  □ 


k=() 


n~l 


(38) 

(39) 

(40) 

(41) 


The  iteration  equations  (25)  and  (26)  define  a  fractional  generalization  of  the  dissipative 
standard  map.  For  3  =  0  this  map  gives  the  Zaslavsky  map  (5)  and  (6)  with 


/'  =  (1  ~c~q)/q 


(42) 


and  fi  =  0. 
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V.  FRACTIONAL  DERIVATIVE  IN  THE  UNKICKED  TERMS  AND  THE  SEC¬ 
OND  FRACTIONAL  DISSIPATIVE  STANDARD  MAP 


In  this  section  we  suggest  a  fractional  generalization  of  the  differential  equation  for  a 
kicked  damped  rotator  with  fractional  derivatives  in  the  uukicked  terms  and  derive  the 
corresponding  discrete  map. 

We  consider  the  fractional  generalization  of  equation  (11)  in  the  form 

OC 

oD?X(t)  -  OoDf.X  (f)  -  ssin(A')  £  6{t  -  «),  (43) 

whore 

q  €  IS,  1  <  a  <  2,  8  sa-1, 

and  o Df  is  the  Riemann  Liouvillo  fractional  derivative  [7  9L  which  is  defined  by  Eq.  (16). 
This  equation  has  fractional  derivatives  in  the  unkicked  terms,  i.e.  on  the  left-hand  side  of 
Eq.  (43).  We  use  the  minus  in  the  left-hand  side  of  Eq.  (43) t  where  q  can  have  a  positive 
or  negative  value. 


Proposition  3.  The  fractional  differential  equation  of  the  kicked  system  (43)  is  equiva¬ 
lent  to  the  discrete  map 


xn+i  =  — — - — —  V  Yk+i  wa{q.  k  -  n  -  1), 

Ya+i  =  **  [>;+esin(Ari) 

where  the  functions  Wa(ayh)  are  defined  by 

W*{atb)  =  ci'"aea^  hva  _  j ,a/>)  -  r(a  -  l,a(fc  +  1)) 

and  F(a,  fc)  is  the  incomplete  Gamma  function  (27). 


(44) 

(45) 

(46) 


Proof.  Let  us  define  an  auxiliary  variable  £(£)  such  that 


c 

0 


Dl 


t  =  *(/•), 


where  ^  D(  °  is  the  Caputo  fractional  derivative  (24).  Using 


0/?”°  =  m  -  m-  (0  <  2  -  a  <  1), 


(47) 


(48) 
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we  obtain 


oD?X  =  Df  o lt~aX  =  Df  o I?~«  $Bf~“t  =  D?m  -f(0))  =  »U, 


-)2  f2  a  C  rV2 


(49) 


and 

a  A*  A  =  A1  o4'flA'  =  D]  olf  "aX  = 

=  A1  off'"1  ?Aa  A  -  Al(€(<)  -  e(0))  =  Ae 

Substitution  of  (49).  (50)  and  (47)  into  Eq.  (43)  gives 

OC 

A2^  -  <iD\i  =  c  sin( £ D2~° 0  ]T  5(t  -  7i).  (1  <  a  <  2). 

n*rO 

The  fractional  equation  (51)  can  be  presented  in  the  Hamiltonian  form 

oc 

Y  —  qY  ~  esin y~\5(/,  -  n),  (1  <  a  <2,  q  e  R). 

n»0 

Using  Eq.  (26)  of  Proposition  2.  we  obtain 


^  «-•  i 

yn+i  =*  +  £sin^— —  VUa(<7,  A;  —  n)^ 


For  (AAni  Vn)>  we  use  equation  (38)  in  the  form 

-  «~»l 

a,.  =  ?  A2;A  =  pTTZY)  E  **+>  k  ~  n)- 

As  a  result,  we  have 


X 


l 


nt-l 


r(«  -  f) 


E^+i  Wa(q,k-n  -  1), 


Jfc-0 


Yn  4 1  —  <Jt  Vn  +csin(X„) 


where  Wa{a,b)  is  defined  in  (46).  This  ends  the  proof.  □ 


If  we  use  the  variables 


(50) 


(51) 


(52) 


(53) 

(54) 


P„  =  fiYn ,  6  = -A  Z  —  —nee? 


11 


Q.  0- 


Om  0 


20 


10 

CX  0 

-10 


K  =  1 2.97.  ft  =  1.9975 


-2^- 


20 

1 

10 

_ 

Q-  0 

^ _ — •»' -  ”  _ _ - — — 

-10 

- -  j 

=  12.93,  a  =  2 

2  X  4 


2  x  4 


FIG.  1:  Structures  of  the  chaotic  attractors  for  different  values  of  K  obtained  after  105  iterations 
(r  —  5.  V.  —  0);  a  =  19975  in  Figs  la-e.  a)  K  =  12  83;  b)  K  —  1 2.87;  c)  K  —  12.93;  d)  zoom  of 
Fig.  lc;  e)  K  =  12.97;  f)  a  —  2,  K  -  12.93. 


then  equations  (44)  and  (45)  give 


t  n 


11  —  p/  ,  .  ,  Pfc+I  11  &  tl  1 ). 

P„.M  =-6Pn-z?8in(Xl)! 


(55) 

(56) 


12 


20 


e 


FIG.  2:  Structures  of  the  fractional  chaotic  attractors  for  K  =  12.93  and  different  values  of  a 
obtained  after  ID5  iterations  (V  =  5,  Q  —■  0). 


These  equations  can  be  considered  as  a  fractional  generalization  of  the  dissipative  standard 
map  equations  (7)  and  (8)  with  Q  =  0.  For  a  =  2,  this  fractional  dissipative  standard  map 
gives  the  dissipative  standard  map  that  is  described  by  Eqs.  (7)  and  (8). 
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CL. 


FIG.  3:  Attracting  trajectories  for  K  ~  12.93  and  small  values  of  a  obtained  after  105  iterations 

(r*5,  n  =  o). 

VI.  NUMERICAL  SIMULATIONS 


Numerical  simulations  were  performed  for  the  second  fractional  dissipative  standard  map 
(Eqs.  (55)  and  (56)).  First  we  used  our  code  to  reproduce  the  results  presented  in  Fig.  1  from 
[19  for  the  structures  of  the  chaotic  attractors  of  the  dissipative  standard  map  at  the  window 
of  the  ballistic  motion  near  K  ^  4t r  ( q  —  —5,  K  =  eexp(q),  and  used  in  [19]  F  is  equal 
to  —q)  for  the  fractional  standard  map  with  a  =  2  and  obtained  a  perfect  agreement  (an 
example  is  given  in  Fig.  If).  As  a  decreases  slightly  from  a  =  2  to  a  =  1.9975,  the  window 
of  the  ballistic  motion  shrinks  and  moves  to  the  higher  values  of  K.  Already  for  a  =  1.9975 
in  Figs,  la-e  the  ballistic  motion  appears  for  K  >  12.86  and  disappears  at  K  ~  12.97.  The 
window  is  completely  closed  at  a  «  1.9969.  The  structures  of  two  symmetric  attractors 
with  disjoint  basins  which  appeal*  within  the  window  (Figs.  lb,c)  is  also  very  different  from 
the  structures  of  the  dying  at  tractors  of  the  dissipative  standard  map  [1.  1 9  The  attractor 
in  Fig.  Id  evolves  from  period  8  trajectory  to  period  4,  period  2,  and,  finally,  period  1 
trajectory  slowly  moving  in  the  direction  of  the  upper  left  corner  with  the  step  of  the  order 
of  10-7. 

When  a  decreases  further,  the  structures  of  the  fractional  chaotic  attractors  evolve  in 
the  manner  presented  in  Fig.  2,  where  one  can  find  one-scroll,  two-scroll,  and  four-scroll 
fractional  chaotic  attractors,  strongly  deviating  from  the  chaotic  attractor  of  the  dissipative 
standard  map  Fig.  If  (see  also  [19]).  The  problem  of  existence  of  multi-scroll  fractional 
chaotic  attractors  was  considered  in  [22  but  for  the  fractional  differential  equations  with 
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the  Caputo  derivatives.  For  values  of  a  near  1  fractional  chaotic  attractor  turns  into  period 
two  and  for  smaller  values  period  one  attracting  trajectories  Fig.  3. 

VII.  CONCLUSION 

The  suggested  discrete  maps  with  memory  are  generalizations  of  the  dissipative  standard 
map.  These  maps  describe  fractional  dynamics  of  complex  physical  systems.  The  suggested 
fractional  dissipative  standard  maps  demonstrate  a  chaotic  behavior  with  a  new  type  of 
attractors.  The  interesting  property  of  these  fractional  maps  is  long-term  memory.  As  a 
result,  a  present  state  of  evolution  depends  on  all  past  states  with  the  weight  functions.  The 
fractional  dissipative  standard  maps  are  equivalent  to  the  correspondent  fractional  kicked 
differential  equations.  Note  that  to  derive  discrete  maps  an  approximation  for  fractional 
derivatives  of  these  equations  is  not  used. 

Computer  simulations  of  the  suggested  discrete  maps  with  memory  prove  that  the  non¬ 
linear  dynamical  systems,  which  are  described  by  the  equations  with  fractional  derivatives, 
exhibit  a  new  type  of  chaotic  motion.  This  type  of  motion  demonstrates  a  fractional  gener¬ 
alization  of  attractors. 

The  special  cases  of  discrete  maps  have  been  studied  to  describe  the  properties  of  the  frac¬ 
tional  chaotic  attractors  of  these  differential  equations.  Under  a  wide  range  of  circumstances 
such  maps  give  rise  to  chaotic  behavior. 
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